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A large class of quintessential inflationary models, recently proposed by Peebles and Vilenkin, 
leads to post-inflationary phases whose effective equation of state is stiffer than radiation. The 
expected gravitational waves logarithmic energy spectra are tilted towards high frequencies and 
characterized by two parameters: the inflationary curvature scale at which the transition to the 
stiff phase occurs and the number of (non conformally coupled) scalar degrees of freedom whose 
decay into fermions triggers the onset of a gravitational reheating of the Universe. Depending upon 
the parameters of the model and upon the different inflationary dynamics (prior to the onset of 
the stiff evolution) the relic gravitons energy density can be much more sizeable than in standard 
inflationary models, for frequencies larger than f Hz. We estimate the required sensitivity for 
detection of the predicted spectral amplitude and show that the allowed region of our parameter 
space leads to a signal smaller (by one f .5 orders of magnitude) than the advanced LIGO sensitivity 
„Cj ■ at a frequency of O.f KHz. The maximal signal, in our context, is expected in the GHz region 

where the energy density of relic gravitons in critical units (i.e. /iq^gw) is of the order of fO 6 , 
roughly eight orders of magnitude larger than in ordinary inflationary models. Smaller detectors (not 
£ — 1 necessarily interferometers) can be relevant for detection purposes in the GHz frequency window. 

£Nj ' We suggest/speculate that future measurements through microwave cavities can offer interesting 

perspectives. 

> 

O . I- FORMULATION OF THE PROBLEM 

o 

The idea that our present Universe could be populated by a sea of stochastically distributed gravitational waves 
(GW) is both experimentally appealing and theoretically plausible. It is appealing since it would offer a natural 
Q\ , cosmological source for the GW detectors which will come in operation during the next decade, like LIGO Q, VIRGO 
0, LISA H and GEO-600 0] R It is also plausible, since nearly all the models trying to describe the first moments 
of the life of the Universe do predict the formation of stochastic gravitational wave backgrounds d|| ■ 

Our knowledge of early the Universe is only indirect. The success of big-bang nucleosynthesis (BBN) offers an 
explanation of the existence of light elements whose abundances are of the the same order in different and distant 
(— > ■ galaxies. BBN hints that when the cosmic plasma was as hot as 0.1 MeV, the Universe was probably dominated by 
radiation jjj. Prior to this moment direct cosmological observation are lacking but one can be reasonably confident 
that the laws of physics probed in particle accelerators still hold. Almost ten years of LEP (Large Electron Postitron 
collider) tested the minimal standard model (MSM) of particle interactions to the precision of the one per thousand 
for center of mass energies of the order of the Z boson resonance. The cosmological implications of the validity of 
the MSM are quite important especially for what concerns the problem of the baryon asymmetry of the Universe and 
of the electroweak phase transition || . In spite of the success of the MSM we have neither direct nor indirect hints 
concerning the evolution of the Universe for temperatures higher than 100 GeV. The causality principle applied to 
the Cosmic Microwave (CMB) phtons seems to demand a moment where different patches of the Universe emitting 
a highly isotropic CMB were brought in causal contact. This is one of the original motivations of the inflationary 
paradigm ||. 

It is not unreasonable to think that in its early stages the Universe passed through different rates of expansion 
deviating (more or less dramatically) from the radiation dominated evolution. It has been correctly pointed out 



Oh 



through the years and in different frameworks |10| that every change in the early history of the Hubble parameter 
leads, inevitably, to the formation of a stochastic gravitational wave spectrum whose frequency behavior can be used in 
order to reconstruct the thermodinamical history of the early Universe. The question which naturally arises concerns 
the strength of the produced gravitational wave background. 



*Electronic address: giovan@cosmos2.phy.tufts.edu 

1 LIGO (Laser Interferometric Gravitational wave Observatory), LISA (Laser Interferometer Space Antenna) 
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If an inflationary phase is suddenly followed by a radiation dominated phase preceding the matter dominated epoch, 
the amplitude of the produced gravitons background can be computed and the result is illustrated in Fig. [|, where 
we report the logarithmic energy spectrum of relic gravitons 



1 dp GW 
p c d In v ' 



(1.1) 



at the present (conformal) time rj as a function of the frequency v ( Pgw is the energy density of the produced 
gravitons and p c is the critical energy density) 0. 

Since the energy spectrum ranges over several orders of magnitude it is useful to plot energy density per logarithmic 
interval of frequency. The spectrum consists of two branches a soft branch ranging between v = 1.1 x 10 -18 ho Hz 
(corresponding to the present horizon) and Udec = 1-65 x 10" 16 (fio ^o) 1 ^ 2 Hz (where Qq is the present fraction of 
critical density in matter and 0.5 < ho < 1 is the indetermination in the experimental value of the Hubble constant). 
For v > i^dcc we have instead the hard branch consisting of high frequency gravitons mainly produced thanks to 
the transition from the inflationary regime to radiation. In the soft branch flowi^, Vo) ~ v~ 2 . In the hard branch 
^gw(^')o) is constant in frequency (or almost constant in the quasi-de Sitter case [see Section VI]). The soft branch 
was computed for the first time in jll) (see also |t2) ) . The hard branch has been computed originally in Jl3| (see also 
@). 
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FIG. 1. We report the graviton spectrum computed in the case of a pure de Sitter phase evolving towards the (present) 
matter dominated phase through an intermediate radiation dominated stage of expansion. The spectrum has a soft branch ( 
10" 18 Hz < v < 10" 16 Hz ) and a hard branch ( 10" 16 Hz < v < 10 9 Hz). In the two black boxes we spot the COBE and the 
millisecon d pu lsar timing bound. Only for illustration purposes we plotted the spectrum for the largest amplitude consistent 
with Eq. (O). 



The COBE (Cosmic Microwave Background Explorer) observations of the first (thirty) multipole moments of the 
temperature fluctuations in the microwave sky imply p5l that the gravitational wave contribution to the Sachs- Wolfe 
integral cannot be larger than the (measured) amount of anisotropy directly detected. The soft branch of the spectrum 
is then constrained and the bound reads 



{l G w{v,Vo)h 2 < 6.9 x 10" 



(1.2) 



for v ~ vq. Moreover, the very small size of the fractional timing error in the arrivals of the millisecond plusar's pulses 
imply that also the hard branch is bounded according to 



(1.3) 



for v ~ 10" 8 Hz corresponding, roughly, to the inverse of the observation time during which the various millisecond 
pulsars have been monitored |H3|. 

The two constraints of Eqs. ( |1.2| ) and (1.3) are reported in Fig. pi at the two relevant frequencies, with black boxes. 
In Fig. |l| we have chosen to normalize the logarithmic energy spectrum to the largest possible amplitude consistent 
with the COBE bound. The COBE and millisecond pulsar constraints are differential since they limit, locally, the 



Notice that in this paper we will denote with In the Neperian logarithm and with log the logarithm in ten basis. 
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logarithmic derivative of the gravitons energy density. There exists also an integral bound coming from standard BBN 
analysis fl7]]l8| and constraining the integrated graviton energy spectrum: 



K ft GW (^7?o)dln^< 0.2 x 1(T 5 , (1.4) 

where y mm corresponds to the (model dependent) ultra-violet cut-off of the spectrum and v n is the frequency cor- 
responding to the horizon scale at nucleosynthesis It should be noted, in fact, that modes re-entering after the 
completion of nucleosynthesis will not incr ease the rate of the Universe expansion at earlier epochs. From Fig. [I] we 



see that also the global bound of Eq. (|L4J) is satisfied and the typical amplitude of the logarithmic energy spectrum 
in critical units for frequencies vj ~ 100 Hz (and larger) cannot exceed 10 -14 . This amplitude has to be compared 
with the LIGO sensitivity to a flat Slcwf"/,*) which could be at most of the order of /iq^gw^/i Vo) = 5 x 10 -11 
after four months of observation with 90% confidence (see third reference in ||). 
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FIG. 2. We illustrate the graviton spectrum produced by a pure de Sitter phase evolving towards a stiff phase. The hard 
branch corresponds to modes leaving the horizon during the de Sitter epoch and re-entering during the stiff phase. The soft 
and the semi-hard branch are made of gravitons re-entered, respectively, during the matter and during the radiation dominated 
phase. As in the previous figure we report the various constraints on the differential spectrum. The frequency at which the 
spike starts developing is model dependent. 



Suppose that the hard branch of the spectrum, reported in Fig. can be split into two further branches, a truly 
hard branch with growing slope and an intermediate semi-hard branch. The situation we are describing is indeed 
reproduced in Fig. || where the semi-hard branch now corresponds to the flat plateau and the hard branch to the 
spike associated with a broader peak. This class of spectra can be obtained in the context of inflationary models 
provided the inflationary phase is followed by a a phase whose effective equation of state is stiffer than radiation. A 
model of this type has been recently investigated in Ref. J2l| by J. Peebles and A. Vilenkin. 

If an inflationary phase is followed by a stiff phase then, as it was showed in 22 1 , one can indeed get a three branch 
spectrum including the usual soft and flat branches but supplemented by a truly hard spike. In general the slope 
of the logarithmic energy spectrum is typically "blue" since it mildly increases with the frequency. More specifically 
the slope depends upon the stiff model and it can be shown 22 that the maximal slope (corresponding to a linear 
increase in J7gw( j/ , Vo)) can be achieved in the case where the sound velocity of the effective matter sources exactly 
equals the speed of light [ p3|p4[| . In Fig. || we illustrate the case of maximal slope in the hard branch corresponding 
to r2G\v(^7%) ~ v\nv. 

Given the flatness of the spectra arising in the case of ordinary inflationary models (see Fig. [l]) the most constraining 
bound comes from large scale observations. In our case the most constraining bounds for the height of the spike and 
for the whole spectrum come from short distance physics and, in particular, from Eq. ( |l.4| ). In order to visually 
motivate the need for an accurate computation of the graviton spectra in the case where an inflationary phase is 
followed by a stiff phase, let us focus our attention on the frequency range where the gravitational wave detectors are 
(or will be) operating. 



3 Notice that the BBN constraint of Eq. ( jhj ) has been derived in the context of the simplest BBN model, namely, assuming 
that no inhomogeneities and/or matter anti-matter domains are present at the onset of nucleosynthesis. In the presence of 
matter-antimatter domains for scales comparable with the neutron diffusion scale [p^|]2c| l this bound might be slightly relaxed. 
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FIG. 3. We draw the hard branch of the relic graviton energy density reported in Fig. |i| with particular attention to 
frequency range where the planned gravitational wave detectors are (or will be) operating. The dashed region marks the 
portion of the spectrum of Fig. |^ which is above the one of Fig. [[]. 

From Fig. || we see that around the LIGO |lj and VIRGO j^j frequency the hard branch of the spectrum has 
a larger amplitude if compared to the case of the spectral amplitude obtained when a pure de Sitter phase evolves 
suddenly towards a radiation dominated epoch. 

In Fig. H we also illustrate (with thick black) boxes the expected sensitivities for interferometric detectors 
(LIGO/VIRGO) and for their advanced versions. In the same figure we also report the expected sensitivities coming 
from the cryogenic, resonant-mass detector EXPLORER, while operating in CERN at a frequency of 923 Hz p5[ . 
EXPLORER provided a bound on r2Gn/<500 which is clearly to high to be of cosmological interest. Nevertheless, by 
cross-correlating the data obtained from bar detectors (EXPLORER, NAUTILUS, AURIGA) it is not unreasonable 
to expect a sensitivity as large as SIgi^ ~ 10~ 4 in the KHz region. Notice finally that the cross-correlation of two 
resonant spherical detectors |2(| might be able to achieve sensitivities as low as Qgw ~ 10~ 7 always in the KHz range. 

Spikes in the stochastic graviton background are not forbidden by observations and are also theoretically plausible 
whenever a stiff phase follows a radiation dominated phase. In this paper, by complementing and extending the 
analysis of |21j and of , we want to study more accurately the spectral properties of the relic gravitons with special 
attention to the structure of the hard peak. We will also be interested in comparing the predictions of the models 
with the foreseen capabilities of the interferometric and resonant detectors. 

The plan of our paper is the following. In Section II we will introduce the basic aspects of the quintessential 
inflationary models. In Section III we will compute the relic graviton energy spectra. In Section IV we will discuss 
the power spectra and the associated spectral densities. In Section V we will compare the obtained spectra with 
the sensitivities of the planned interferometric and non-interferometric detectors. In Section VI we will analyze the 
impact of the slow-rolling corrections on the structure of the hard peak. Section VII contains our concluding remarks. 
For sake of completeness we made the choice of reporting in the Appendix some relevant derivations of the formulas 
used in obtaining our results. 



II. QUINTESSENTIAL INFLATIONARY MODELS 

Recently J. Peebles and A. Vilenkin [^lj presented a model where the idea of a post-inflationary phase stiffer than 
radiation is dynamically realized. One of the motivations of the scenario is related to a recent set of observations 
which seem to suggest that f^o (the present density parameter in baryonic plus dark matter) should be significantly 
smaller than one and probably of the order of 0.3. If the Universe is flat, the relation between luminosity and red-shift 



observed for Type la supernovae 27 seem to suggest that the missing energy should be stored in a fluid with negative 
pressure. The missing energy stored in this fluid should be of the order of 10 -47 GeV 4 , too small if compared with 
the cosmological constant arising from electroweak spontaneous symmetry bereaking (which would contribute with 
(250 GeV) 4 ). The idea is that this effective cosmological constant could come from a scalar field 4> (the quintessence 
p8| field) whose potential is unbounded from below |2j|. According to Peebles and Vilenkin, could be identified 
with the inflaton and, as a result of this identification, the effective potential of <f> will inflate for (f> < and it will 
be unbounded from below for </> > acting, today, as an effective (time dependent) comological term. A possible 
potential leading to the mentionad dynamics could be 

AM 8 

y{4>) = H4> 4 + M 4 ), for0<O, and V(cf>) = ^ - — 4 , for > 0. (2.1) 
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where, if we want the present energy density in <j> to be comparable with (but less then) the total (present) energy 
density, we have to require M ~ 10 6 GeV. The scenario we are describing can be implemented with any other 
inflationary potential (for cj> < 0) and the example of a chaotic potential is only illustrative. Our considerations will 
be largely independent on the specific potential used and we will comment, when needed, aboout possible differences 
induced by the specific type of potential. 

Let us consider the evolution equations of an inflationary Universe driven by a single field (f> in a conformally flat 
metric 

ds 2 = a 2 (77) (dr/ 2 - dx 2 ). (2.2) 
Using the conformal time 77 the coupled system describing the evolution of the scale factor and of <fi is 

Mpi 2 = + a 2 V{<jy) 

MpH' = - (d> n - a 2 V(4>) 

((>" + 2H(f>' + a 2 — =0, (2.3) 



where l P = Mp 1 = ^8ttG/3. Since the scalar field potential is unbounded from below, after a phase of slow-rolling 
the inflaton evolves towards a phase whe re t he kinetic energy of the inflaton dominates. For instance one can bear in 



mind the form of V((j)) reported in Eq. (2.1). The background enters then a stiff phase where the energy density of 



the inflaton and the scale factor evolve as 



H 2 M 2 (—), o(»j)~.E where, H x = — ~ VXM P . (2.4) 



771 a\r)x 

When an inflationary phase is followed by a stiff phase a lot of hard gravitons will be generated. At the same time 
the energy density of the background sources will decay as a~ e whereas the energy density of the short wave-lengths 
gravitons will decay as a~ 4 . Te Universe will soon be dominated by hard gravitons whose non-thermal spectrum p3] 
would be unacceptable since gravitons cannot thermalize below the Planck scale. A solution to this potential difficulty 
came from L. Ford ]3(| who noted that in the limit of nearly conformal coupling also scalar degrees of freedom (possibly 
coupled to fermions) are amplified. If N s minimally coupled scalar field are present they can reheat the Universe with 
a thermal distribution since their energy spectra, amplified because of the transition from the inflationary to the stiff 
phase, can thermalize thanks to non-graviational (i.e. gauge) interactions which get to local thermal equilibrium well 
below the Planck energy scale. It can be also shown that the same discussion can be carried on in the case where the 
scalar degrees of freedom are simply non conformally coupled |3l| . 

Suppose indeed that during the inflationary phase various scalar, tensor and vector degrees of freedom were present. 
Unless one adopts some rather contrived points of view we have to accept that, in Einsteinian theories of gravity, 
the only massless degrees of freedom to be amplified by a direct coupling to the background geometry are tensor 
fluctuations of the metric and non conformally coupled scalar fields since the evolutions equations of chiral fermions, 



gravitinos |32| and gauge fields p3j are invariant under a Weyl rescaling of the metric tensor in a conformally flat 



background geometry as the one specified in Eq. (2.2). Of course, if the theory is not of Einstein Hilbert type this 
statement might be different. 

The evolution equation of a non conformally coupled scalar field in a conformally flat FRW background reads 

ip" + 2Hip' - 6£ [H' + H 2 ] - V 2 ip = 0. (2.5) 

By defining the corresponding proper amplitude x = a ^ we get that the previous equation can be written, in Fourier 
space, as 

X 'k+[k 2 -V(v)]xk = 0, V( V ) = (l-60 — , (2.6) 

a 

where we see that the case of exact conformal coupling is recovered for £ — > 1/6 whereas the case of minimal coupling 
occurs for £ — > 0. A lot of work has been done in the past in order to compute the energy density of the quanta of the 
field ip, excited as a result of the background geometry evolution in the early stages of the life of the Universe |m| . 
One can try to do the calculation either exactly (but only for rather specific forms of the effective potential of the 
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Schroedinger-like equation ( J2.6| )), or approximately by identifying |£ — 1/6| as the small parameter in the perturbative 
expansion. In this limit one can show |p0| , j34| that the energy density of the created quanta can be expressed as 

MV) = J l^( fc )| 2 fc3 with P( k ) =2kJ_ e'^VWdri. (2.7) 

We can notice that in most of the examples we are interested in V{rj) — > ?/~ 2 for 77 — > ±00. As an example we can 
consider V{rj) — q/(rj 2 + J7 2 ). Then by evaluating f3(k) with contour integration in the complex plane we can estimate 
that the energy density will be 

^ 128 a 4 128 u V ' 

where 7ii = 17/771, and Hi = Tii/a is the Hubble parameter in cosmic time. Suppose now that during the inflationary 
phase there are N s (minimally coupled and massless) scalar degrees of freedom ipi. Because of the minimal coupling 
to the geometry these scalar degrees of freedom will clearly be excited since their evolution equation is not invariant 
under conformal rescaling of the metric tensor. The produced quanta associated with each ipi can be computed by 
specifying (for each field mode) the initial vacuum state deep in the de Sitter epoch and by ensuring a sufficiently 
smooth transition between the de Sitter and the stiff phase. We will perform a similar calculation for the case of 
GW in the next Section. Here we only report the main result which was originally obtained in |3(if| in the case of 
quasi-conformal coupling (i.e. |£ — 1/6| < 1) and subsequently generalized to the case of generic £ in pl[ : 

M«?) = E^)~«Hf(-) . where R~RiN 8 . (2.9) 

Ri is the contribution of each massless scalar degree of freedom to the energy density of the amplified fluctuations 
and it is of the order of 10 -2 . The appearance of H\ in the final expression of the energy of the created quanta 
can be simply understood since the typical spectra obtained in the transition from a de Sitter phase to a stiff phase 
are increasing in frequency p2| ] and, therefore, the maximal contribution to the energy density will come from the 
ultra-violet branch of the spectrum. 

The creation of massless quanta of the fields ipi triggers an interesting possibility of gravitational reheating. Since 
p m red-shifts faster than p^, we have that there will be a moment, r\ r where the two energy densities will be of the 
same order. In the context of the present model this moment defines the onset of the radiation dominated phase. We 
can compute this moment by requiring that p m (j]r) ~ P4>{Vr)- The result is that 



^)^V!j|^~io-Vii, (2.10) 



a r J V Mp 

where we used the fact that in order to be compatible with the COBE observation \f\ ~ 10~ 7 |55|. In view of our 
application to GW it is interesting to compute the typical (present) frequency at which the transition to radiation 
occurs. By red-shifting the curvature scale at ?7 r (i.e. H r = H(rj r ) = H^ai/a,.) 3 ) from r) r up to now we obtain 

*fo,)=3.58iZ* (^) (^) 1/3 mHz, (2.11) 
-LU yth 

where gth is the number of spin degrees of freedom contributing to the thermal entropy after matter thermalization. 
Amusingly enough this frequency is of the same order of the typical frequency of operation of LISA. It is also interesting 
to compute the present value of the frequency v\ , with the result that 

i^(Tjb) = 358 R~ * (^£) 1/3 GHz. (2.12) 
5th 

The thermalization of the created quanta of the fields ipi occurs quite rapidly, and its specific time is fixed by the 
moment at which the interaction rate becomes comparable with the Hubble expansion rate during the stiff phase. The 
typical energy of the created quanta is of the order of e ~ Hi(ai/a). The particle density is of the order of n ~ Re 3 . 
Assuming that the created quanta interact through the exchange of gauge bosons, then the typical interaction cross 
section will be of the order of a ~ a 2 /e 2 . Thus, imposing that at thermalization 71(77*) (7(77*) ~ H(rj±) wc get that 
a±/ai ~ a _1 i? -1 / 2 , with a ~ 10 _1 -10 -2 . The typical temperature associated with the transition from the stiff to the 
radiation dominated phase can be computed and it turns out to be 



G 



Tr = (J±L)R^M P ~ 10 3 N^ 4 GeV. 



(2.13) 



If we do not fine-tune H\ to be much smaller than 10 -7 in Planck units and if we take into account that N s has to 
be typically large in order to be compatible with standard BBN (see also Section III) we have to conclude that T r is 
typically a bit larger than 1 TeV. 



III. GRAVITONS ENERGY SPECTRA 



We can characterize a generic graviton background in terms of three related (and equally important) physical 
observables. We can compute the (present) spectral energy density in critical units Hcwt^ ^o), but, for experimental 
applications, two other quantities can be defined, namely the power spectrum ( which will be denoted with 5h(u, rj )) 
and the spectral density Su{v, rjo). flcwi^, vo) and 8h(y, Vo) are dimension-less whereas the spectral density is measured 
in seconds. In Appendix A we give the precise mathematical definitions of these observables. 

The continuity of the scale factors and of their first derivatives implies that the evolution of our model can be 
expressed as 



<H(ri) = 

a s (v) = 
a r (v) = 

a m (r]) = 



m 
v 



2r] + 3r?i 



m 

77 + 3?7i 4- T) r 
^Vl( 2 Vr + 3?7i) ' 

(v + Vdcc + 6?7i + 2r\ r ) 2 



V < -m, 

- m < r) < rj r , 
Vdcc <V <r)r, 

VO < V < Vdcc, 



(3.1) 



where 770 and fydec are, respectively, the present time and the decoupling time, whereas 771 and r\ r have been defined 
in the previous Section. 

The graviton field operators can be decomposed as 



1 



(27T)* 



ik-: 



(3.2) 



where fi(k,rj) — if)(k,rf)a(k). This decomposition holds for each polarization. In order to compute the energy density 
of the graviton background we have to solve the evolution of the mode function 



V/' + [fc 2 



u = o, 



(3.3) 



in each of the four temporal regions defined by Eq. (3.1). Notice that a" /a has a bell-like shape and it goes 
asymptotically as rf 2 in each phase of the background evolution. Thus tp will oscillate for krj ^> 1 but it will be 
parametrically amplified in the opposite lim it ( i.e. krj < 1). At krj ~ 1 the giv en mode will hit the potential barrier 
represented by \a"/a\. The solution of Eq. (3.3) in the background of Eq. (3.1) is: 





v Ik 






V < -r)i, 


l()s(k,7l) = 




(y) + sA_ 




- ?7i < Tj < Tj r , 


ip r (k,rj) = 


l ^\B + {k)e~ lz + £L 
V2k l 


(k)e% 




Vr < V < Vdcc, 


^m{k,7j) = 


-j=V^[q*c+(k)H{? 


(w) + qc- 


-(A)ff«H], 


Vdcc < V < VO, 



where, 



P = 



^ 2 ^g 1 4 



(3.4) 



(3.5) 
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(1 2) 

guarantee that the large argument limit of the Hankel functions H^i ' is exactly the one required by the quantum 
mechanical normalization The arguments of tp are, respectively, 



ki-j, y = k(i] + -r/i), z = krj, w = k(r/ + r/ dcc + 6r]i + 2rj r ), 



(3.6) 



and the six mixing coefficients (A±(k), B±(k), c±(k)) can be fixed by the six conditions obtained matching ip and 
ip' in r\ = —Tfi, rj = rj r and 77 = ?7dcc- The results of this calculation are reported in Appendix C. For a gene ric 
amplification coefficient (3-{ui) the spectral energy density in the relic graviton background is given by Eq. (A. 22) 



dpow 
d\nuj 



7 n(uj), n{ui) = |/3_(w)| 2 , uj 



k 
a 



2-nv 



(3.7) 



since, as it is well known and discussed in Appendix A, the square modulus of the mixing coefficient can be interpreted 
as the mean number of gravitons at a given frequency. Notice that uj is the physical wave-number. The relic gr avito n 
energy spectrum (in c ritical units) in each o f the three branc hes is simply obtained by inserting A-{uj) (i.e. Eq. ( ]C.2p , 
B-(ui) (i.e. Eq. (C.7) and C-(uj) (i.e. Eq. ( C.1C )) into Eq. ( A. 22 ). The final result can be expressed as 



n GW (oj,j] ) = fl 7 (r]o) e A 



UJ r 



hi' 



n GW (oj,7] ) = n 7 (»?o) 4 e A ln 

7T 

!Jgic(w,i)o) = ^yivo) Yg e A 



In 2 ^ 



UJ r < UJ < UJi, 



^doc < UJ < LO r , 



UJO < UJ < Wdcc, 



with 



^dcc(ryo) = 1-69 x f(T 16 x (fioM/io) 1/2 Hz ' and 
where fio^o) is the fraction of critical density in matter. Notice that 

'5doc\l/3 „ 81 



1.1 x 1(T 18 ho Hz, 



(3.8) 



(3.9) 



gth 



n 7 (77o) 



Pc(vo) 



go 



32 7T 3 

rpi 
1 



30 HlM 2 p 



2.6 x 1CT 5 h^ 2 , 



(3.10) 



where go = 2 and To = 2.73 K. ^(770) is the fraction of critical energy density in the form of radiation at the present 
observation time. Notice that the dependence upon the number of relativistic degrees of freedom occurs since, unlike 
gravitons, matter thermalizes and then the ratio between the critical energy density and the energy density stored in 
the relic graviton background is only approximately constant in the radiation dominated phase. 

The local (differential) bounds on the energy spectrum can be easily satisfied. Indeed by taking Hi/Mp < 10 -7 the 
spectrum satisfies the COBE bound of Eq. (1.2) and also the pulsar bound of Eq. (|1.3|). The indirect nucleosynthesis 



bound applies to the integrated spectrum and since in our case the spectral energy density increases sharply in the 
hard branch we have to conclude that the height of the peak cannot be too large. In order to prevent the Universe 
from expanding too fast at nucleosynthesis we have to demand 



7 

d\nuj£l GW (uj,r) D ) < ^(N v - 3) 



Pc(Vn) 



(3.11) 



Since the maximal number of massless neutrinos permitted in the context of the homogeneous and isotropic BBN 
scenario is bounded to be N u < 3.4, we have that in our context the nucleosynthesis bound becomes 



3 

~N~s 



5n_ 

3th 



1/3 



< 0.07, 



(3.12) 



Notice that we kept the Hankel indices v and fj, generic. In the case of a pure de Sitter phase we would have v 
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where the factor of 3 counts the two polarizations of the gravitons but also the quanta associated with the inflaton. 
The number of relativisti c deg rees of freedom after matter thermalization is given, in the MSM by gth = 106.75, 
whereas g n = 10.75. Eq. ( |3.12 ) implies that the number of (minimally coupled) scalar de gree s of freedom will have 
to exceed 20 as it can occur, for instance, in the minimal supersymmetric standard model pij . 

An increase in N s does not only decreases the height of the peak, but it can also make the peak structure narrower. 
This happens simply because by increasing N s , v r oc grows and v\ oc N s gets pushed towards more infra-red 
values. Given the limited range of variation of R this effect is quite mild We illustrate the variation of R on the energy 
spectrum in Fig. || A decrease in the inflationary curvature scale at the end of inflation does not affect the peak since 
the maximal amplified frequency does not depend on Hi/Mp but it does only depends on R. 




FIG. 4. In the left plot we report the energy spectrum as a function of the physical wave number for a fixed value of A ( 
~ 1(F 7 ) but for different values of R. In particular we show the cases R — 10 (dot-dashed line), the case R — 1 (full thick 
line) and the case R = 0.6 (full thin line, almost invisible). We see that by increasing R the height of the peak gets smaller 
and also its width shrinks. In the right plot we report the graviton energy spectrum for a fixed value of R (which we take of 
order one) but for different values of A and, more specifically, A = 10" 14 (full thick line), A = 10" 15 (full thin line), A = 10" 16 
(dot-dashed line). We took fio = 0.3, <?dec = 3.36 and gth = 106.75 as fiducial set of parameters. Notice that in spite of the 
fact that Qgw(u, go) does depend on the specific value of ho which appears, in our notations, in fi 7 (01,170). We choose to plot 
/ij(!gw(w, 770) which does not depend upon ho- In this way the the amplitude of the energy spectrum is independent on ho and 
the only quantities depending on ho are uio and lo^cc- In this plot took ho — 0.6. 



In quintessential inflationary models the energy density of relic gravitons can be much larger, at high frequencies, 
than in the case of ordinary inflationary models where the energy spectrum is still flat in the hard branch. The 
location of the peak is rather surprising. In fact it depends (very weakly, as we said) on the number of minimally 
coupled scalar fields but it does not depend upon the final curvature scale at the end of inflation. Thus the peak is 
firmly localized around 100 GHz and it cannot move of one order of magnitude. This behavior has to be contrasted 
in ordinary inflationary models where the maximal frequency of the spectrum is determined by 10 11 y/Hi/Mp Hz. 
So by lowering H\, then, the maximal frequency decreases. In the case we are discussing Hi only appears in the 
expression of v r . Therefore, by decreasing H\/Mp (see Fig. || right plot) vi does not move but v r gets comparatively 
smaller reducing the frequency range of the semi-hard branch. 



IV. GRAVITONS POWER SPECTRA AND ASSOCIATED SPECTRAL DENSITIES 



The relic graviton spectrum can be characterized not only in terms of the energy density but also in terms of 
the power spectrum, namely in terms of the Fourier transform of the two-points correlation function of the graviton 
field operators. Of course the energy density and the power spectrum can be precisely related. On the basis of the 



derivation reported in Appendix A (see Eq. (A. 25)) we can connect the power spectrum to the energy density of the 
relic graviton background 

n GW (w,»jb) = A(^) 2 I^(^%))I 2 , \s h (k, m )\ 2 = fc 3 (IMMo)l 2 + IMMo)l 2 ), (4.1) 

07T Ho' 

where h®(k, if) = il>®{k,rf)/a{rf) and h§{k,rf) — ip®(k, rf)/a(rf) are t he F ourier amplitudes of the graviton field operators 
associated with the two (independent) polarizations [see Eqs. ( jA.Sj )] and Hq is the present value of the Hubble 
parameter. The three branches of our power spectrum turn out to be: 
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6 h (u>,r)o) = 6.5 X 10 J r ; ; In — , w < w < w doc 4.2 

8 V iio (t7o ) V w / W/ 

where £> = 3.07 x 10~ 30 /io- Our results are illustrated in Fig. |[ As we can see the power spectrum of the hard 
branch evolves typically as uj^ 1 / 2 . Our power spectrum declines slower than in ordinary inflationary models where 
the high frequency tail evolves typically as a-> _1 . This behavior occurs, in our case, for frequencies Wdcc < uj < 0J r . 
By varying R in the allowed range the power spectrum is only slightly affected. We stress this point in Fig. |5| where 
different power spectra are reported for various values of R at fixed lambda. 
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FIG. 5. At the left we plot the hard branch of the power spectrum for a fixed value of A (which we take of order 10 _1 ) 
but for different values of R and, more specifically, R = 1 (full thick line), R = 10 (full thin line), i? = 0.1 (dot-dashed line). 
At the right we plot the power spectrum for a fixed value of R (which we take to be 10) but for different values of A and, more 
specifically, A = 10 -14 (full thick line), A = 10 -15 (dot-dashed line), A = 10~ 16 (full thin line). In this and in the previous plot 
we took ho = 0.6 and fio = 0.3 as fiducial values. 



A similar effect can be observed if we keep R fixed and we let A free to change. Again the impact of the variation 
of A affects the power spectrum less than the energy density. The reason for this behavior is due to the fact that a 
variation in 5h{ui,rjo) (at a given frequency) boils down to a quadratic variation in Q,qw{<jJ, rjo). 

In comparing the produced graviton spectrum with the experimental sensitivities of the various detectors it turns 
out to be useful to translate the physical information contained into the energy density into another quantity, the 
spectral amplitude (often called also spectral density) whose relation with the energy density has been derived in 
Appendix A: 

Hz 2 

S h {v, Vo) = 8 x 10~ 37 n GW {v, no) hi (4.3) 



From Eq. (4.3) S) l {u 1 tjq) turns out to be quite a small fraction of a second. For this reason sometimes the experimental- 



1/2 

ists express their bounds in terms of SJ (y, rio). The maximal signal expected from a stochastic graviton background 
is limited (from above) by the nucleosynthesis bound. Let us then demand that the peak of the graviton spectrum 
does not exceed the nucleosynthesis bound. Then, to reach a level of sensitivity comparable with /iqSIqw^ Vo)^i 10~ 6 
will imply that 

V^(^)<3 x 10- 26 0*¥\ V! Hz- 1 / 2 . (4.4) 

The physical relevance of Sh{v,T]o) is related to the way we hope to observe in the near future stochastic GW 
backgrounds. In order to detect a gravitational wave background we in an optimal way |3(| we need at least two 
detectors (two bars, two interferometers, one bar and one interferometer...). Suppose then that we have two detectors 
and suppose that the output of a the detectors is given by Si — hi + A/i where i = 1,2 refers to each single detector 
; hi is the gravitational fluctuation to be detected and Mi is the noise associated with each detector measurement. 
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Now, if the noises of the two detectors are not correlated, then, the ensemble average of the Fourier components of 
the noises is stochastic, namely 



W{v)M i {i/)) = -8(v-i/)8 ii Sj*{v), 



(4.5) 



where Stf{v)is the spectral density of the noise. T he no ise level of the detector can then be estimated bt Sj\f(v). 
As we discussed in Appendix A (see Eqs. ( A..26 )— ( A.27 )) it is also possible to characterize the signal with the same 
technique and, then, we will have 



(hi(v)h*(v')) = -8ij5(v - v')Sh{v)- 



(4.6) 



where Sh{v) = Sh{—v) is the spectral density of the signal and it is related to Hgw by Eq. (4.3). Very roughly, if a 
signal is registered by a detector this will mean that Sh(v*,r]o) > SVo„)' namely, the spectral density of the signal 
will be larger, at a given frequency than the spectral density of the noise associated with the detector pair. In 
order to confront our signal with the available sensitivities we need to compute the spectral density Sh- Taking into 
account the numerical factors, we get 
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Sh(u,r)o) 
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Hz" 



LOQ < U> < Wdec 



(4.7) 



where C ~ 2.12 10~ 13 B 2 . 

Not only the spectral amplitude of the theoretical signal depends upon the frequency, but also spectral amplitude 
of the noise does depend upon the frequency. It is not only important if Sh{v) > Sjj{v) at a particular frequency but 
it is also crucial to take into account, for detection strategies, the spectral behavior of the signal versus the spectral 
behavior of the noise in the frequency range explored by the detectors. The spectral density in the hard branch is 
illustrated in Fig. [| from to = 1 Hz until u> = 1 KHz. We remind that this range of wave-numbers is the one relevant 
for the forthcoming interferometric detectors. 
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FIG. 6. We plot the spectral density for A = 1CP 1 in the frequency range relevant for the forthcoming interferometric data 
on direct GW detection. In the full thick line we report the case R = 10, in the full thin line the case R — 1. and in the 
dot-dashed line the case R = 0.6. Again we took ho = 0.6 and fio = 0.3 as fiducial values. 



In Fig. |6J we see that the spectral density of our signal is mainly concentrated in the blank region between the thick 
line (corresponding to R = 10) and the full thin line (corresponding to R = 1). For to ~ 0.1 kHz, Sh ~ 10 -52 -10 -53 
sec. For oj ~ 0.01 kHz, Sh ~ 10 -50 -10~ 51 sec. This observation shows that, within the frequency range of the 
interferometers our theoretical signal can be larger or smaller depending upon the frequency. 
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V. DETECTABILITY OF THE QUINTESSENTIAL GRAVITON SPECTRA 



There are, at the moment various interferometric detectors under construction. They include the two LIGO 
detectors being built by a joint Caltech/MIT collaboration, the VIRGO detector (near Pisa, Italy) j| the GEO-600 
(Hannover, Germany) Q and the TAMA-300 (near Tokyo, Japan) S|. The noise spectral densities of these detectors, 
defined in a frequency range going from 1 Hz to 10 4 Hz, decline usually quite rapidly from 1 to 100 Hz, they have 
a minimum (around 100 Hz) corresponding to the maximal sensitivity and then they rise again with a more gentle 
slope until, approximately, 1-10 kHz. 

Noise Power Spectra 

(for the major interferometers) 
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FIG. 7. The predicetd noise power spectra for various interferometers. This Figure is adapted from j37]]. P(f) is what we 
called Sjif and it is the quantity which should be compared with Sh- In this plot / = v. 



As we discussed in the previous Section, in order to have some hopes of detection we have to demand that the 
theoretical spectral density of the signal is larger than the spectral density of the noise. Let us then try to compare this 
value with the expected sensitivity of the interferometers possibly available in the near future. The power spectrum 
of our signal is too small to be seen by TAMA-300, GEO-600, and VIRGO. By correlating different detectors the 
sensitivity can increase also by a large factor |j^-E]J. However, the published results on the foreseen sensitivities at 
2/j arc far too large to be relevant for our background j26|. By comparing Fig. |] with Fig. |?] we can argue that only 
the advanced LIGO detectors are closer to our predicted spectral density and that our signal in generally smaller than 
the advanced LIGO sensitivity. The two (identical) LIGO detectors are under construction in Handford (Washington) 
and in Livingston (Lousiana). After various years of operation the detectors will be continuously upgraded reaching, 
hopefully, the so-called advanced level of sensitivity. 

Let us estimate the strength of our background for a frequency of the order of 0.1 kHz -1 kHz. Let us assume that 
the energy density of the stochastic background is the maximal compatible with the nucleosynthesis indications. The 
graviton energy density (in critical units) at a frequency vi ~ 0.1-1 kHz is then 

H G wh,?)o) h 2 Q = 2.29 10~ 15 N~ 3/4 [-19.7 + 0.25 lnA s ] 2 , Vl = 0.1 kHz, (5.1) 
t! GW (^,r) ) h 2 a = 2.29 10~ 14 N~ 3/4 [-17.4 + 0.25 lnN s ] 2 , v T = 1 kHz. (5.2) 

Suppose then that we correlate the two LIGO detectors for a period r = 4 months. Then, the signal to noise ratio 
(squared) can be expressed as f|l| 

S\ 2 _ 9H 4 _ /^..^M^wf^ 



The function "f(u) is called the overlap function. It takes into account the difference in location and orientation 
of the two detectors. It has been computed for the various pairs of LIGO-WA, LIGO-LA, VIRGO and G EO- 600 
detectors [fHf . For detectors very close and parallel, 7^) = 1. Basically, 7^) cuts off the integrand of Eq. (5^3) at 



a frequency 27rz/ of the order of the inverse separation between the two detectors. For the two LIGO detectors, this 

(1 2) 

cutoff is around 60 Hz. SjJ are the noise spectral densities of the two LIGO detectors and since the two detectors 
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are supposed to be identical we will have that Sj^^) = Sffl \v). In order to detect a stochastic background with 90% 
confidence we have to demand S/N>1.65. No w in order to estimate the signal to noise ratio we need to estimate 
numerically the integral appearing in Eq. (5.3). We know the theoretical spectrum since we just estimated it. The 
noise spectral densities of the LIGO detectors are not of public availability so that we cannot perform numerically this 
integral. In the case of a flat energy spectrum we have that the minimum Ogw detectable in r = 4 months is given, 
with 90 % confidence, by OgwC^i, Vo) = 5x 10~ 6 /i^ 2 (for the initial LIGO detectors) and by f^ow(^, %) = 5x lO -11 /^ 2 
(for the advanced LIGO detectors). For a correct comparison we should not confront Eqs. (5.1) and (5.2) w ith the 
sensitivities of the LIGO detectors to flat energy spectra b ut r ather with the sensitivities obtained from Eq. (5.3) in 
the case of our specific energy spectrum reported in Eqs. (3.8). 



(2)/ 



Let us then compare, for illustration Eqs. (5.1 



and (5.2) with the sensitivity to a flat energy spectrum e ven if thi s 
is not completely correct. The idea is to discuss, at fixed frequencies, the maximal signal provided by Eqs. (5.1)-(pM^) 
for different values of R. 
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FIG. 8. With th e tw o full thin lines we illustrate log h^Q,Gw{vi , T)o) as a function of logiV s for frequencies vj = 0.1 kHz 
according to Eqs. ( |5.l[ )-(^^). The full thick lines represents the sensitivity of the LIGO detectors and of the advanced 
LIGO detectors to an energy density with flat frequency spectrum. From top to bottom the thick lines refer, respectively, to 
/ijSJgw = 5 x 10 -6 and to /ijSJgw = 5 x 10 -1 . In order to be detected, the theoretical signal has to be above the thick line(s). 
We see that the region between the two thin lines does not overlap with the sensitivity of the upgraded LIGO detectors by, 
roughly, 1.5 orders of magnitude. 



This comparison is illustrated in Fig. |g. For the allowed range of variation of N s our signal lies always below (of 
roughly 1.5 orders of magnitude) the predicted sensitivity for the detection, by the advanced LIGO detectors, of an 
energy density with flat slope. The main uncertainty in this analysis is however the spectral behavior of the sensitivity 
for a spectrum which, unlike the one used for comparison, is not flat. It might be quite interesting to perform this 
calculation in order to see which is the precise sensitivity of the LIGO detectors to a spectral energy density as large 
as 10 -12 and rising as (y/v r ) In [v/v{) in a frequency range 1 Hz— 1 kHz. 

If we move from the kHz region to the GHz region the signal gets much larger than in the ordinary inflationary 
models. Moreover, the energy spectrum exhibits a quite broad peak whose typical amplitude can be as large as (but 
smaller then) 0.5 x 10~ 5 . The spike (corresponding to the maximum of the peak) is located at a frequency of the 
order of 350 iT 1 / 4 GHz. 

On top of the interferometric detectors there are other types of detectors which are operating at the moment or 
might be operating in the future. The available upper limits on the strength of a stochastic background of relic 
gravitons come from the (cryogenic) resonant-mass detectors. In particular EXPLORER, while operating at CERN, 
provided an upper (direct) bound on the graviton energy density. The frequency vb relevant for this type of detectors 
is among 907 and 923Hz. By analyzing the 1991-1994 data of the EXPLORER antenna, the Rome group of G. 
Pizzella got a bound on the strength of the gravitons energy density of the order of /iq^gw^300. This bound, as it 
is, is not crucial since we know that the graviton energy density should be much lower, however we stress that this 
is direct bound from an operating device. Moreover, by correlating two bar detectors with the same features (like 
EXPLORER, NAUTILUS or AURIGA Q) it is not excluded that a very interesting sensitivity of /i^GW ~ 1CU 3 will 



NAUTILUS (in Frascati, Rome) and AURIGA [|42| (in Legnaro, Padova, Italy) are both resonant mass detectors. 
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be reached. 

The particular spectral shape of the signal coming from quintessential inflation seems to point towards the use 
of electromagnetic detectors and, in particular of microwave cavities. A typical signature of the background we are 
discussing in the present paper is that the peak frequency (which almost saturates the nucleosynthesis bound) occurs 
for frequencies of the order of V\ = 3.58 x 10 11 i? -1 / 4 Hz. We can say that, very roughly, the size of the GW detectors 
is not only determined by construction requirement but also by the typical frequency range of the spectrum we ought 
to explore. In this sense the large separation between the two LIGO detectors is connected with the fact that the 
explored frequency range is of the order of 100 Hz. Thus, if we deal with frequencies which are of the order of the 
GHz we can expect small detectors to be, theoretically, a viable option. 

Microwave cavities can be used as GW detectors in the GHz frequency range [^,[45| . These detectors consist of an 
electromagnetic resonator, with two levels whose frequencies v s and v a are both much larger than the frequency j^gw 
of the gravitational wave to be detected. In the case of [Q the two levels are achieved by coupling two resonators one 
symmetric in the electric fields and the other antisymmetric. Indeed, in the case of cylindrical microwave cavities there 
are different normal oscillations of the electric fields. In the relevant mode for the experimental apparatus is the 
TEon according to the terminology usually employed in electrodynamics in order to identify normal modes of a cavity 
corresponding to different boundary conditions |44j . There were published results reporting the construction of such 
a detector 46 1. In this case j/qw = 10 GHz and Av = v s — v a ~ 1 MHz. In this experiment a sensitivity of fractional 
deformations 8x/x of the order of 10~ 17 was observed using an integration time At ~ 10 3 sec. The sensitivity to 
fractional deformations can be connected to the sensitivity for the observation of a monochromatic gravitational wave 
of frequency vgw- Following we can learn that the sensitivity to fractional deformations is a function of P a and 
P s ( the powers stored in the symmetric and antisymmetric levels) , Q (the quality factor of the cavity [Q and which 
gives rate of dissipation of the power stored in the cavity). If we would assume, as in Q ~ 10 9 , P a ~ 10 -21 Watt 
P s = 2.4 x 10" 2 Watt , we would get Sh ~ lO" 17 . 

There, are at the moment, no operating prototypes of these detectors and so it is difficult to evaluate their sensitivity. 
The example we quoted jh3| refers to 1978. We think that possible improvements in the Q factors can be envisaged 
(we see quoted values of the order of 10 12 which would definitely represent a step forward for the sensitivity). In 
spite of the fact that improvements can be foreseen we can notice immediately that, per haps, to look in the highest 
possible frequency range of our model is not the best thing to do. In fact from Eq. (^ij) we can argue that in order 
to detect a signal of the order of /iqSIqw ~ 10~ 6 at a frequency of 1 GHz, we would need a sensitivity of the order of 
Sh ~ 10~ 30 . Moreover as stressed in jiTj the thermal noise should be properly taken into account in the analysis of 
the outcome of these microwave detectors. Indeed, as noticed from the very beginning fig ], the thermal noise is one 
of the fundamental source of limitation of the sensitivity. An interesting strategy could be to decrease the operating 
frequency range of the device by going at frequencies of the order of 1 MHz. Based on the considerations of ]43j ] 
we can say that by taking high quality resonators the foreseen sensitivity can be as large as /iq^gw ~ 10~ 4 . This 
sensitivity, though still above our signal, would be quite promising. 



VI. GRAVITON SPECTRA FOR QUASI-DE SITTER PHASES 

During an inflationary phase the evolution is not exactly of de Sitter type. We want to understand how the slope 
of the energy density of the relic graviton background will be modified by the slow-rolling corrections for frequencies 
accessible to the forthcoming interferometers. In general the deviations from a de Sitter stage can be induced either 
because of the specific inflationary model or because of the slow-rolling corrections whose strength can be described 
in terms of the so-called slow-rolling parameters 

« = -_< 1, [3 = -^-< I, (6.1) 

H(j> 

where H = (In a)' is the Hubble parameter in cosmic time, cj> is the inflaton and the dot denotes derivation with 
respect to cosmic time. In the slow rolling approximation the inflaton evolution is dominated by the scalar field 
potential according to the (approximate) equations 

ay 

3#</> + — ~ 0, MlH 2 ~ V. (6.2) 
acp 

In this approximation H is not exactly constant but it slowly decreases leading to what we called quasi de Sitter 
phase. The evolution equation of the mode function can be written, by using the definition of a as 

V/'+[fc 2 -2a 2 tf 2 (l-|)]V; = 0. (6.3) 
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As we can see, the time-dependent frequency appearing in the mode function contains two contribution: the first 
one ( i.e. 2a 2 H 2 = 2/i] 2 ) is the term coming from a pure de Sitter phase, the second one, proportional to a, is the 
correction. 

In short the logic is the following. The quasi de Sitter phase modifies (through and a-dependent correction) the 
index of the Hankcl functions whose precise value ( equal t o 1 .5 in the pure de Sitter phase) gets slightly smaller 



than 1.5. From the sign of the correction appearing in Eq. (3J5) we can argue that the quasi de Sitter nature of the 
inflationary phase will lead to a decrease in the slope of the energy spectrum. The question is how much the slope 
of the hard branch will be affected, or, more precisely, how much smaller than one will it be the slope of the hard 
branch of the graviton energy density. 

The answer to this question will of course depend upon the specific inflationary model since the size of the slow- 
rolling corrections can vary from one inflationary potential to the other. Notice that our concern is different from the 
one usually present [[48|-|5~i| in the context of ordinary inflationary models where the slow-rolling corrections are taken 
into account in the soft branch of the spectra (namely at sufficiently large scales). Indeed, for flat spectra the most 
significant bounds come from the infra-red, whereas, in our case the most significant bounds are in the ultraviolet and 
we have to understand how the scales which re-enter in the stiff phase are affected by the quasi-de Sitter nature of the 
inflationary phase. Of course, we could also discuss the slow-rolling coprrections to the soft branch of our spectrum, 
but, they are, comparatively, less relevant (for the structure of the spike) than the correction to the hard branch. 

As usual, the slow-rolling corrections not only affect the mode function evolution but also the definition of conformal 
time itself, namely we will have 

/dt f da 1 f da , . 

-a=\^B=-^\ a ^B- (6 - 4) 

By using the fact that MpH = — (3/2)0 2 and the fact that f3 < 1 we can connect directly a to the slope of the 
potential 



H _ Mp fd\iiV 
H ~ ~6~ 



(6.5) 



Since we are interested in the lowest order slow-rolling correction we will assume that a and (3 are constants. This 
is a simplification which will not affect (numerically ) the slope of the spectrum. In the case of inflationary models 
with chaotic and exponential potential it can be shown that the "running" of a with </> (and, therefore, with the 
wavenumber k) will affect the spectral slopes with a term which is of the order of 1/N 2 where N is the number of 



Eq. (I6J) 



inflationary c-folds 51 53] . If a is constant then we have that 77 1 = — (aH)/(l + a) which leads, once inserted in 



v 2 - i 3 
ip" + \k 2 5-4-1 = 0, with v = - + a, (6.6) 

T] Z 2 

where the expression of v holds for a < 1 . 

Let us examine now the quasi-de Sitter nature of different inflationary scenarios. An inflationary potential of chaotic 
form 

V{4>) = ^, (6.7) 

will lead to a = (Mp/Q)(n 2 /4> 2 ). Let us take the value of a the one corresponding to modes crossing the horizon 
around 20-25 e-folds before the end of inflation. The reason for this choice is very simple. We want to understand 
how is modified the slope of the energy spectrum for scales larger than (and of the order of) the ones probed by the 
interferometers. A simple calculation shows that, for instance, the LIGO/VIRGO scale crosses the horizon roughly 
21 e-folds before the end of inflation. The end of inflation occurs when a{4> cn d) = 1- Thus, at the end of inflation 



'/'end = (Mp/6)n 2 . By solving consistently Eqs. (6.2) with the potential ( |6.7| ) we get easily that the number of 
inflationary e-folds at a given value of 4> is given by 



da 

From this last equation we can determine easily 021 and then, by inserting <fi 2 nd , a turns out to be 

n 

n + 84' 



(6.9) 
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which is of the order of 0.02 for n — 2 , of the order of 0.04 for n = 4 and so on. Another example could be the one 
of an exponential potential. Using the definition of a from Eq. (|6.5| ) we have that for an exponential potential of the 
form 

U(0)-exp[iAl ; (6.10) 



P" 



will lead to a = l/q 2 - 

With these results we can easily compute the corrections to the hard branch of the spectrum. The spectral energy 
density will be 

^Gw(w,f?o) = n 7 (i7o)/(a, X,R)( — V 2 "ln 2 (— ), with 

f(a,X,R) = A ( a +|) 2 2 2 « |r(§ + atfX^Br**^) 1 ' 3 . (6.11) 

The power spectrum and the associated spectral density can be computed from using the techniques already discussed 
in the previous Sections. Some authors call —2a the corrections to the spectral index. Eq. ( |6.11 ) tells us that the 



corrections arising during a quasi-de Sitter phase always go in the direction of making the maximal slope of the hard 
branch slightly smaller than one by a factor 0.08 (for instance in the case the case of 4> 4 potential). In the case of 
exponential potential the magnitude of the correction to the slope is controlled by q 2 . Again we see that for reasonable 
values of q 2 like q 2 = 5, 10, ... || the correction are again small and, at most, of the order of few percents. 



VII. CONCLUDING REMARKS 



In this paper we computed the relic graviton spectra in quintessential inflationary models. We showed that the 
energy spectra possess a hard branch which increases in frequency. Large energy densities in relic gravitons can be 
expected in this class of models. The spike of the energy spectrum can be as large as /iqIIqw — 10~ 6 at a typical 
(present) frequency v\ = 358 x i?" 1 / 4 GHz. The spectral amplitude at the interferometers frequencies is just below the 
value visible by the upgraded LIGO detectors. Since a large amount of energy density can be stored around the GHz 
the use of small electromagnetic detectors for the detection of such a background seems more plausible. In particular 
microwave cavities should be considered as a possible candidate. Our investigation also hints that the sensitivities of 
the advanced LIGO detectors to energy spectra increasing, in frequency, as u\nv should be precisely computed by 
convolving our spectra with the noise power spectra of the detectors. 
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APPENDIX A: RELIC GRAVITONS CORRELATION FUNCTIONS 



In this Appendix we concentrate some of the more technical derivations concerning the spectral properties of the 
relic gravitons which can be characterized through their energy density in critical units or through the two point 
correlation function of the amplified tensor fluctuations of the metric. We want to give explicit derivations of these 
quantities and of their relations. 

Let us start from the canonical action describing relic gravitons in a curved metric of Friedmann- Robertson- Walker 
type. GW, being pure tensor modes of the geometry, only couple to the curvature but not to the matter sources. 
Their effective action can be obtained by perturbing the Einstcin-Hilbert action 



to second order in the amplitude of the tensor fluctuations of the metric 

9p,v (x, rf) -> g MV {rf) + 8gfJ {x , rf) , 5g liv = h liv , 



0. 



(A.l) 



(A.2) 



Pure tensor fluctuations of the metric (with two physical polarizations) can be constructed by using a symmetric 
three-tensor hij satisfying the constraints h\ — 0, Vj/ij = (where Vj is the covariant derivative with respect to 
the three-dimensional background metric). The action perturbed to second order in the amplitude of the tensor 
fluctuations becomes 



5 (2) S (T) 



(A.3) 



We can construct two independent scalars corresponding to the two polarization states of the gravitational waves: 

hij(x,ri) = h 9 {x,r))q^'{x) + h®(x,r))q^> (x), (A.4) 

where 



9(i) (x) 



\{x)e{{x) - ei(x)el(x)}, qf 2) {x) = -[e\{x)e{{x) + 4{x)e\{x)}. 



(A.5) 



In these coordinates e±, €2 and e*3 form an orthonormal set of vectors. In the previous decomposition we assumed the 
wave to propagate along e^. So the perturbed action can also be written, in units of 61 p — 1, as 



6 (2) S (T) = 1 / ' <p Xy pjj[d a h e d l 3he& xf> + d a h®d{3h®g a ^- 

The associated energy-momentum tensor can be written as [p5|.p(i[ 



d^hsdvh® + dph^dvh® - g [ g a Pd a h@dph® + g a Pd a h®dph 



The canonically normalized graviton field can be written as Hij = ahij and its corresponding action reads 



= \J 



d xdrj 



(A.6) 



(A.7) 



(A. 



With this form of the action we can promote the classical canonical normal modes to field operators obeying the 
standard (equal time) commutation relations. The Fourier expansion of the field operators will then look 



1 



(2tt) 3 / 2 
1 

(2tt) 3 / 2 



d 3 k 
d 3 k 



Ae(fc,77)e* 5 + /i| B (fc,77)e-* ? 



(A.9) 



where jj,^{k,rj) = 77)0® (fc) and jj,®(k,r]) = ip®(k,rj)a®(k). The creation and annihilation operators follow the 

usual commutation relations namely 
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[a x (k),aUp)]=S x ^ 3 \k-p). 



(A.10) 



From the Hciscnbcrg equations of motion for the field operators follow the evolution equations for the Fourier ampli- 
tudes 



1>" + [k 2 



(A.ll) 



for each of the two physical polarizations. In most of the cosmological applications we will be dealing with \a" /a\ 
has a bell-like shape going to zero , for large conformal times, as i]~ 2 . Eq. ( A.llj ) can be be solved in two significant 
limits. The first one is for k 2 <C \a" /a\ corresponding to wavelength which are outside of the horizon (under the bell) 



ip{k,rj) = A(k)a(r)) + B(k)a(rj) 



(A.12) 



where A(k) and B(k) are integration constants. For large momenta (i.e. k 2 3> \a"/a\) the general solution of Eq. 
(|A.11|) is 



i>(k, V ) = -±= [0+ {k)e-^ + /?_ [k)e^] 
\/2k 



(A.13) 



where (3-{k) and (3+{k) are complex numbers and where the quantum mechanical normalization l/\/2k has been 
chosen. Therefore, if we start at r) — > — oo with a positive frequency mode the evolution for k 2 <C |a"/ a l wm m i x 
the the positive frequency with the negative one and eventually we will end up, for r\ — ► +00 with a superposition of 
positive and negative frequencies according to Eq. (A.13) 

The two point correlation function of the graviton field operators can then be computed in Fourier space 

(O|/*A(M)/4(a»0|O) = \MKri)\ 2 h^\k-p). (A.14) 

Field modes with different wave-numbers are then statistically independent as a consequence of the graviton emission 
from the vacuum. Some authors refer to Eq. (A.14) as to the stochasticity condition which express the main statistical 
property of the graviton background. Within our quantum mechanical formalism we can also compute the two-point 
correlation function between field operators, namely 



Z(r) = (0\( llJ (x)^(x + r)\0) = ^ 



dk sin kr 
k kr 



\Sh(k,v)f 



where 



|<MM)l 2 = ^(l^©(M)l 9 + Mfe(M)l a ), 



(A.15) 



(A.16) 



is the power spectrum of the relic gravitons background summed over the physical polarizations. 

Sometimes, to facilitate the transition to the quantities often used by the experimentalists, it is convenient to denote 
the gravitational wave amplitude, for each polarization, with 



h(x,r)) = 



1 



(2tt) 3 / 2 



d 3 k h(k,r)) e 



(A.17) 



where the reality condition implies h*(k) — h(—k). The statistical independence of waves with different wave vectors 
implies, in this formalism, that 



(h x (k, V )h*(p,r,)) = \h x (k,r,)\ 2 S X(T 6^(k-p), h x (k, V ) = ^M. 



(A.18) 



Eq. (A.18) can be viewed as the classical limit of Eq. (A.14) where the quantum mechanical expectation values are 
replaced by ensemble averages. Notice that this concept can be made quite rigorous by taking the average of the 
graviton field over the coherent state basis. This procedure would be analogous to what is normally done in quantum 
optics in the context of the optical equivalence theorem |54| . 

The energy density (and pressure) of the relic graviton background can be easily obtained by taking the average of 
the energy momentum tensor. The energy density is 
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Pgw(v) = (0|T 00 |0), 

which can also be written as 

Pgw(v) = (O^OelToolOeO^) = 



(A.19) 



1 



167r 3 a 2 



d A k 



K(k, V )\ 2 + K(k, V )\ 2 + k 2 [\h e (k,r,)\ 2 + |MM)I 



(A.20) 



If we insert the as ympt otic expression of the mode function of each polarization for large (positive) conformal times 
we get, from Eq. ( A..20 ), 



PGwyV) = / — w = , w = fe/a, 

7 w n 2 



from which we can define the logarithmic energy spectrum as 

dp GW (uj,r]) uj 4 2 
jj = -2 P-M , 



(A.21) 



(A.22) 



(where we used the fact that \f3+(u>)\ 2 — \/3-(oj)\ 2 = 1) By taking into account that the energy density of n(uS) gravitons 
in the proper momentum interval duj is given by 



dpaw — 2ujn(uj) 



d 3 LU 

(2^)3' 



(A.23) 



|/3_(fc)| 2 can be interpreted as the mean number of produced gravitons in a given frequency interval. 

In order to compute the relation between Q,Qw{u),ri) and 5h(u},r)) — 5^(u>, rf)/a{rf) we can use Eq. ( A. 17 ). The 
energy density in gravity waves is given by 



PGW{V) = ^-2\ h <g> 



h: 



/2 



(A.24) 



where the brackets indicate now ensemble average of the Fourier amplitudes. Using now Eq. (A. 17) together with 
the stochasticity condition we obtain (after having divided by the critical energy density) 



n GW (u,r]Q) = (-^-) 2 |<5 h (w,?7o)| 2 - 



(A.25) 



where the rjo is the present conformal time which appears since we want to refer everything to the present value of 
the Hubble constant. In comparing the signal coming from a particular model with the experimental sensitivities it 
turns out to be widely used the spectral density 5^(^,770). In order to define it let us assume that each polarization 
of the gravity wave amplitude can be written as 

i7rut , (A.26) 



(A.27) 



h{v) = J dt hit) e~ 2 " v \ 
where t is the cosmic time variable. Then the spectral amplitude can be defined as 

{h{v)h*{v')) = \6{v-v')S h {v). 



By inserting this expression into Eq. (A.24) and by using the definition (A.27) we obtain 

47T 2 

This last equation implies that 



OUr. 



S h {v,r) ) = 8 x 10- 67 n G w{v,r]o)K 



(A.28) 



(A.29) 



from which is clear that Sh is measured in seconds. 



19 



APPENDIX B: RELIC GRAVITONS SPECTRA IN THE ORDINARY INFLATIONARY CASE 



Consider the model of an Universe evolving from a de Sitter stage of expansion to a matter-dominated phase passing 
through a radiation dominated phase. By requiring the continuity of the scale factors and of their first derivatives at 
the transition points we have that (1(77), in the three different temporal regions, can be represented as 



<H(i]) = 
a r {rj) = 
a m (r)) = 



Tjl 
?7 + 27/1 



for 
for 



(v + m + ^m) 2 

4t7i (772 + 2?7i) 



v < -m 

- Tji < T] < r/2 
for 77 > 772 



(B.l) 



where —771 and r/2 mark, respectively, the onset of the radiation-dominated phase and the decoupling time. 

In order to compute the graviton spectra in this model we h ave to estimate the amplifi cation of the graviton mode 
function induced by the background evolution reported in Eq. (B.l). The solution of Eq. (A. 11) in the three temporal 
regions is given by 



Vv(y) 



P 



2k 

[B+e 



B-e*»], 



VZ[q*c+H}?\z)- 



where 



kr], y = k{f] + 2t7i), z = k(rj + r] 2 + iffr). 



(B.2) 



(B.3) 



(1 2) 

and Hv A are the Hankel functions 



)8| of first and second kind. In the case of the background given by Eq. (B.l) 
the Bessel indeces v and /z are both equal to 3/2 but we like to keep them general in light of the applications reported 
in the present paper. 

Notice that in Eq. (|B.2|) we included 



P = 



= Jl p -ifa+2u) g = 1 /![ fi <}(l+2M) 



(B.4) 



ensuring that the large time limit of the mode function is the one dictated by quantum mechanics (i.e. e l±kr, /\/2k) 
without any extra phase or extra (co nsta nt) coefficient. 

The constants appearing in Eq. (B.2) are fixed by the quantum mechanical normalization imposed during the 
de Sitter phase, and then, in order to compute the amplification we have to match the various expressions of the 
mode functions (and of their first derivatives) in r\ = —r\\ and in r\ = 772- By doing this we get the expression of the 
amplification coefficients. By firstly looking at modes which went outside of the horizon during the de Sitter phase 
and re-entered during the radiation dominated phase (i.e. modes 27r/?7i < k < 2^/772) we have that B± are given by 



^=(" + 5)]±«V=5T^2i(-» 1 )} 



(B.5) 



We are interested in the _B_ which measures the amount of mixing between positive and negative frequency modes 
and whose square modulus can be interpreted as the mean number of produced gravitons. Moreover, we want to 
evaluate 73_ in the small argument limit (i.e. x\ -C 1) which is the one physically relevant. The result is 



B.(k) = 



I» 



-if (1+2*/) 



I kr)i 



(B.6) 



where T(v) is the Eulcr Gamma function. In the case of a pure de Sitter background we get, indeed, 



B-(k) = - 



-\krjx 



(B.7) 



20 



In order to compute the spectrum of the relic gravitons crossing the horizon during the matter dominated epoch 
we have to consider a further transition by matching ip r (y) and ip m (z) (and their first derivatives) in 77 = 772. The 
result is that 



c- = -7.(5+ - B_)^Hl 2 \z 2 )q, (B. 



where we reported only the dominant part in the X\ -C 1 and z 2 *C 1 limit. By performing the limit explicitly we get 

\c-(k)\ = r( ^J (l/) 2"+^ - ^)MNN x r- * N 2 r+i (B.9) 



In the case of the model of Eq. (B.l) we have that the energy density of the created gravitons in critical units at 
the present observation time 

^0) = -^, (B.10) 
p c amv 



can be computed by inserting into Eq. (A. 22) the expressions of the mixing coefficients of Eqs. (B.6)-(B.8), with the 
result that 

u \ 2 / \ 2 



This spectrum is reported in Fig. The so/t branch and the Aard branch defined in the Introduction do correspond, 
respectively, to the two frequency ranges v$ < v < v d cc and v& cc < v < v\. 

APPENDIX C: TRANSITION FROM INFLATION TO STIFF PHASE: ACCURATE MIXING 

COEFFICIENTS 



In order to determine the six matching coefficients appearing in Eq. ( |3.4D we have to match the mode function ip 
and its first derivative in r\ = —771, 77 = r] r and 77 = ?7dcc- The results of this calculation are reported in the present 
section. Consider first the amplification leading to the hard branch of the spectrum. The mixing coefficients are given 
by 

A + (A0 = ^e-i-{ff«(^)[a^ (C.l) 

The small argument limit of A- (directly relevant for the estimate of the graviton spectrum in the hard branch) can 
be easily obtained and it turns out to be 

A-(k) ~ — 2 v -i e -^ {2v+1) T(v)x^\Yixi. (C.2) 

In order to derive the previous and the following expressions it is useful to bear in mind that the Wronskian of the 
Hankcl functions is given by 

W[H^(r),H^(r)] = H ( ^ +1 {r)H^(r) <{r)H {2) +1 {r) = -ii (C.3) 

for a generic argument r and for a generic Bessel index a. In this calculation it is also useful to recall that the small 
argument limit of the Hankcl functions is 
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where the minus (plus) refers to the Hankel function of first (second) kind. Notice that this formula is only valid for 
a^O. In this last case we have indeed that 



JZ< 1,2) (r)~lT-lnr. 



(C.5) 



Let us now consider the following transition namely the one leading to the semi-hard branch of the spectrum. The 
exact expression of the mixing coefficients is, in this case, 



B±(k) = ^{y/y;[ s *A + (k)H^\y r )+ S A-(k)H^\y r )] ± i[s* A + {k){^H^ {y r ) - ^.H?\y r )) 
+sA_ (k) (^=H^ (y r ) (»r ))] } • 



The small argument limit of B_ leads to 



B_(fc) = -e 



which becomes, in the case v = §, 



-^=2-1^)1^1^ \k m \- In [f], 



(C.6) 



(C.7) 



(C.8) 



Finally, let us compute the amplification coefficients in the case of the soft branch of the spectrum. In this case we 
have that the mixing coefficients are 

C _(fc) = l{(B + (fc)e-^« +B_(fc)e < ^)[( At + J) -§L#f (*dec) 

Z Z ^/ Zdcc 

-V^«*^2i(2dec)] +*(B+(fc)e-^ - J B_(fc)e^)Vi^ g *^ 2 )(^ d ec)}. 
c_(fc) 

-V^^i^dcc)] +i(B+(fe)e- i! ' d '» -B-e^-^v^gff^^dcc)}. 
In the small argument limit we have 

c_(fe) =-^2^-i e -^+ 2 ^r( M )rH|fc^r^|A : r 7l rnfcr; dcc |^+hn [^-], 

47T rj r 



l -{{B + e-^ +B_{k)e i ^)[{n+ \) -L= H^(z dec ) 

Z Z a/ Zdcc 



which becomes, in the case fi = 3/2, 

c-(fc) = 



9 7T 



^^e-^l^r^lAfMecl-^l-'In [^]. 



(C.9) 



(CIO) 



(C.11) 



The derivations we just reported are the basis for the results presented in Section III. 
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